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We minimize the dissipation rate of an electrical network under a global constraint on the sum 
of powers of the conductances. We construct the explicit scaling relation between currents and 
conductances, and show equivalence to a a previous model [J. R. Banavar et al Phys. Rev. Lett. 
84, 004745 (2000)] optimizing a power-law cost function in an abstract network. We show the 
currents derive from a potential, and the scaling of the conductances depends only locally on the 
currents. A numerical study reveals that the transition in the topology of the optimal network 
corresponds to a discontinuity in the slope of the power dissipation. 
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The optimal distribution of valuables such as electric- 
ity or telephone signals has been a subject of much study 
since Westinghouse and Edison's War of the Currents in 
the late 19*h century ^f|; more recently, natural systems 
such as river networks and vascular systems have been 
fruitfully interpreted in this light Q, 0, • Hence formal 
models of optimal transport networks have attracted at- 
tention over many years However, different studies 
use different definitions of network and optimize differ- 
ent functionals. For example, Durand [y,l3 considers hy- 
draulic networks whose currents derive from a potential, 
explicitly analogous to electrical networks; the networks 
are embedded in an ambient space, and he studies the 
optimal geometry and the relation between the local ge- 
ometryand local topology. On the other side, Banavar 
et al. 3] propose a more abstract model where the graph 
is not assumed to be embedded in a target space, nor are 
the currents through the nodes explicitly constrained to 
derive from a potential. This allows them to furnish a 
strictproof that the topology of the optimized flow pat- 
tern [3j depends on the convexity of their cost function, 
but makes a direct physical interpretation of the model 
more elusive. In the following, we shall introduce a third 
model of an optimal transport network from whom both 
of these previous models can be derived, so all formula- 
tions are, in fact, equivalent. 

Consider an electrical transport network on a graph com- 
posed of nodes k interconnected by links (k, I). There is a 
given current source if. at each node and the total current 
input must add to zero: ^2 k ik — 0. There are variable 
currents Iki flowing through the links; the sum of all cur- 
rents impinging on a given node k must equal the given 
current sources: ik = Y]i Iu (Kirchhoff's current law). 
We associate a resistor Rki > to each link (k, I) and 
decompose its value as Rki = (dkiKki) -1 , where dki > 
is a given weight and the conductances km are variable; 
considering fcfcl CIS HI conductivity per unit length, d k i can 
be thought of as the length of the link. The dissipation 



rate J is then a function of the currents Iki through the 
links and the conductances km- 
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We shall minimize this dissipation rate J over the cur- 
rents Iki and the conductances Km with the local con- 
straint given by Kirchhoff's current law, and a supple- 
mentary global constraint that the sum over the conduc- 
tances raised to a given power 7 > is kept constant: 



tf 7 = E*« 



One may interpret this constant as an amount of re- 
sources we have at our disposal to build the network. 


Since we allow km and Iki to vary independently, the 
currents are not explicitly constrained to derive from a 
potential at the nodes Uk and Kirchhoff's voltage law 
(the sum of the potential differences on a loop vanishes) 
need not apply. 

Using a Lagrange multiplier A, we define the function 




FIG. 1: Sketch of a loop a indicating the direction of the cur- 
rents. Every perturbation of the Iki satisfying the constraints 
can be written as a weighted sum of such loops. 
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FIG. 2: Examples of the conductivity distributions. Results of the relaxation algorithm with different initial conditions: (a) 
7 = 2.0 and (b) 7 = 0.5. (c) Result with optimized topology with 7 = 0.5. the arrow indicated the localized inlet, the remaining 
nodes are outlets with constant 
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The necessary conditions for a minima of J with constant 
K are then: 



= 0, 
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Let us first consider the derivatives with respect to Jjy. 
Let {/fcz}, minimize J. Adding a circular current 

JT Q on a loop a to the currents (fig. 0) does not violate 
the constraints. We (re)define the directions of the cur- 
rents Jfez on the loop to be parallel to the loop current 
X*. Then 

dz 



loop Q 

Thus Kirchhoff's voltage law holds at the minimum of 
J, so the currents though the links derive from potential 
differences between the nodes: Iki = Rki{Ui — Ufs)- Note 
that this is not the case for every arbitrary current dis- 
tribution. For instance, if all currents on the loop in fig. 
^are positive Iki > 0, then there exists no set of Rki > 
to fulfill this relation. 

Let us now consider the derivatives of S with respect to 
Kki (eq. EJ). With the constraint of a constant K, we 
obtain an explicit scaling relation between the currents 
and the conductivity in the minimal configuration: 
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We can now write the total dissipation (eq. ^) in terms 
of the currents alone as 
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Since for 7 > 0, the function x 1+ ~ is monotonically in- 
creasing, the original minimization problem is reduced to 
the minimization of 



C({/w» = E( i *«/ d «) 



By setting 
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and rescaling the weights as Wki = d kl 1+7 , the quantity 



to be 



zed is 
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which is exactly the model used by Banavar et al.^. 
They give a strict proof that for T < 1, the resulting 
structure may not have any loop, and each spanning tree 
is a local minimum. For T > 1, there are in general loops 
and a unique minimum. Due to the correspondence be- 
tween 7 and r, this result must apply also to our original 
model where 7 < 1 (7 > 1) corresponds to a T < 1 

(r>i). 

On the other hand, the correspondence between the dif- 
ferent models allows an important conclusion about the 
model of Banavar et al. Since in both formulations, the 
minimum is obtained by the same set of currents, and 
since in our model these currents must derive from po- 
tential differences between the nodes, this must be true 
for the minimum of the Banavar et al. model, too. We 
can furthermore write down directly the values of the 
corresponding resistors as 



Rkl = (dklKkl) 
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with an arbitrary positive constant A. Rki thus scales 
explicitly with the local currents for V 5^ 2. 
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Since positive 7 corresponds to < T < 2, the equiv- 
alence of the two models is restricted to this parameter 
range. T > 2 corresponds to values 7 < — 1, for which 
our model collapses into infinitely many degenerate min- 
ima. The relations correspond instead to a saddle node 
of J: a minimum with respect to the Im and a maximum 
with respect to the km- Nevertheless direct inspection 
shows that the current flow in the Banavar et al. model 
is potential with the set of resistors given by eq. 1101 even 
for r > 2. [13 

In order to get a deeper insight into the transition at 
7 = 1, we search numerically for the minimal dissipa- 
tion configuration of an example network, a triangular 
network of conductivities with a hexagonal border, with 
equal weights df-i = 1. The total number of nodes N no d es 
scales roughly as the square of the linear dimension of 
the network, given by the diameter of the graph N^i a - 
Except for those on the border, each node is linked by 
conductivities to six neighboring nodes. 
We place a current source at a corner of the hexagon (io), 
the remaining (N noc ies — 1) nodes present homogeneous 
distributed sinks; each node absorbs i/- = —io/{N no 4 es — 

1)- 

As an order parameter, we will consider the normal- 
ized dissipation rate J mini J homo-, where J homo is the to- 
tal dissipation with a constant conductivity distribution 
Kki = const., and J m in is the dissipation for the opti- 
mized distribution of the conductivities. Note that 3 homo 
corresponds also to 7 — > 00. 

The previous discussion allows us to simplify the mini- 
mization problem enormously: using the scaling relation 
between km and Im, one can restrict the search of the 
minimum to the space of the currents or the space of 
conductivities. Furthermore, we can use the fact that 
the optimized current distribution derives from a poten- 
tial Uk to construct a simple relaxation algorithm. Start- 
ing with a random distribution of km, we calculate first 
the values of the potential at the nodes by solving the 
system of linear equations ik = J2i Rki{Uk — U{), then 
the currents through the links Im are determined. We 
use these currents to determine a first approximation of 
the optimal conductivities on the basis of the scaling re- 
lation. Then, the currents are recalculated with this set 
of conductivities, and the scaling relation is reused for 
the next approximation. These steps are repeated until 
the the values have converged. We check by perturbing 
the solution that it actually is a minimum of the dissipa- 
tion, which was always the case. 

For all 7 > 1, independently of the initial conditions, the 
same conductivity distribution is obtained, which con- 
forms to the analytical result of [3( : there exists a unique 
minimum which is therefore global. 

Furthermore, the distribution of km is "smooth" , varying 
only on a "macroscopic scale", as show in Fig.^a)). No 
formation of any particular structure occurs. However, 
the conductivity distribution is not isotropic. We can 
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FIG. 3: (a) The normalized minimum dissipation rate 
Jmin I ' J homo as a function of 7 for a network with Ndia = 15 
(462 links) and a network with Ndia = 31 (2070 links, in 
red). Note the discontinuity of the slopes at 7 = 1. (b) A 
detailed view of the crossover at 7 = 1 for Ndia ~ 15. Cross 
symbols show data points obtained by optimizing a tree topol- 
ogy, circles show the output of the relaxation algorithm. The 
continuous lines indicate the actual minimum. 



interpret the conductivity distribution as a discrete ap- 
proximation of a continuous, macroscopic conductivity 
tensor (see also 6]). The smooth aspect of the distri- 
bution is conserved while approaching 7 — > 1 while the 
local anisotropy increases, while the values of all km re- 
main finite, even if they get very small. For 7 = 1.5 and 
Ndia = 15, the conductivity distribution spreads already 
over eight decades and becomes still broader as 7 — * 1 + , 
in which limit the number of iteration steps diverges as 
the minima becomes less and less steep. 
7 = 1 presents a marginal case. The results of the sim- 
ulation suggest that the minimum is highly degenerate, 
i.e., there are a large number of conductivity distribu- 
tions yielding the same minimal dissipation. 
For 7 < 1, the output of the relaxation algorithm is quali- 
tatively different (fig-|2tb))- The currents are canalized in 
a hierarchical manner: a large number of conductivities 
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rapidly converge to zero and thus vanish transforming the 
topology from a highly redundant network to a spanning 
tree. This, too, is predicted by the analytical results 
In contrast to 7 > 1, the conductivity distribution can 
not be interpreted as a discrete approximation of a con- 
ductivity tensor: for Ndi a — ► 00, the structure becomes 
fractal. 

For different initial conditions, the relaxation algorithm 
yields trees with different topologies: each local minima 
in the high-dimensional and continuous space of con- 
ductivities {nu} correspond to a different tree topology. 
Given a tree topology, the currents through the links are 
given directly by the topology and do not depend on 
the values of the km, and so using the scaling relation, 
one may directly write down the dissipation rate for a 
given tree. For 7 < I, we do thus not need to apply 
the relaxation algorithm, but we should search for the 
global minima in the (exponentially large) space of tree 
topologies using a Monte-Carlo algorithm. This regime 
has been widely explored in the context of river networks 
0,0,0], mainly for a set of parameters that corresponds, 
in our case, to 7 = 0.5. An example of a resulting mini- 
mal dissipation tree structure is given in fig- EI C ) • Note 
also, that the scaling relations can be seen as some kind 
of erosion model: the more currents flows trough a link, 
the better the link conducts 

The qualitative transition is reflected also quantitatively 
in the value of the minimal dissipation (fig[Sfa)). The 
points for 7 > I were obtained with the relaxation 
algorithm, the points 7 < 1 by optimizing the tree 
topologies with a Monte-Carlo algorithm. For 7 — > 00, 

J rain I J homo 

— ¥ 1 by definition; for 7 —> 0, J m in/Jh omo r 
0, because the vanishing km allow the the remaining 
Kkl —> 00. 

Figure Bib) shows the behavior of minimal dissipation 
rate close to 7 = 1. For 7 smaller than one, the relax- 
ation method spends a long time only to furnish a local 
minimum, while the Monte-Carlo algorithm searching for 
the optimal tree topologies gives lower dissipation values. 
The different values corresponding to different realization 
indicate that the employed Monte-Carlo method does not 
find the exact global minima. For 7 > 1, the relaxation 
algorithm gives the lower J m i n because the global min- 
ima does not have a tree topology. 

While the curve is continuous, the crossover at 7 = 1 
shows a clear change in the slope of Jmin(?()- One could 
interpret this behavior as a second order phase transi- 
tion. (The change in slope is of course preserved in the 
function C(Im) used by |3j-) 

As an intriguing practical application of these models, 
one may for instance cite the venation of plant leaves. 
Experimental evidence jjj shows that the water transport 
through the veins derives from a pressure gradient. The 
venation pattern however shows a enormous redundancy 



of loops 0, 0, E|- On the basis of some examples, it 
has been proposed 0, 0] that the loops are actually 
meaningful to optimize the water transport in the leaf. 
The results presented in this paper however shows that 
this is not the case: optimization cither leads to a tree 
topology, or to no structure at all. If the venation pattern 
is really based a optimization principe, it cannot simply 
be optimization of a steady state water transport, even if 
Murray's law seems to hold at the nodes of the venation 
0. 
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